Mon. Not. R. Astron. Soc. 00Q.fTHT7lf201 1) Printed 2 March 2013 (MN KTeX style file v2.2) 



Predicting the Fate of Binary Red Giants Using the Observed 
Sequence E Star Population: Binary Planetary Nebula Nuclei and 
Post-RGB Stars 



J. D. Nie 1 ' 2 *, P. R. Wood 2 * and C. P. Nicholls 3 * 

1 Department of Astronomy, Beijing Normal University, Beijing, 100875, China 

^Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, Weston Creek ACT 2611, Australia 
3 Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA 92093, USA 



Accepted; Received; in original form 



ABSTRACT 

Sequence E variables are close binary red giants that show ellipsoidal light variations. 
They are likely the immediate precursors of planetary nebulae (PNe) with close binary cen- 
tral, stars as well as other binary post-AGB and binary post-RGB stars. We have made a Monte 
Carlo simulation to determine the fraction of red giant binaries that go through a common en- 
velope (CE) event leading to the production of a close binary system or a merged star. The 
novel aspect of this simulation is that we use the observed frequency of sequence E binaries in 
the LMC to normalize our calculations. This normalization allows us to produce predictions 
that are relatively independent of model assumptions. In our standard model, and assuming 
that the relative numbers of PNe of various types are proportional to their birthrates, we find 
that in the LMC today the fraction of PNe with close binary central stars is 7-9%, the fraction 
of PNe with intermediate period binary central stars having separations capable of influencing 
the nebula shape (orbital periods less than 500 years) is 23-27%, the fraction of PNe contain- 
ing wide binaries that are unable to influence the nebula shape (orbital period greater than 500 
years) is 46-55%, the fraction of PNe derived from single stars is 3-19% and 5-6% of PNe 
are produced by previously merged stars. We also predict that the birthrate of post-RGB stars 
is ~4% of the total PN birthrate, equivalent to ~50% of the production rate of PNe with close 
binary central stars. These post-RGB stars most likely appear initially as luminous low-mass 
helium white dwarf binaries. The average lifetime of sequence E ellipsoidal variability with 
amplitude more than 0.02 magnitudes is predicted to be ~0.95 Myr. We use our model and the 
observed number of red giant stars in the top one magnitude of the RGB in the LMC to pre- 
dict the number o f PNe in the LMC. We predict 548 PNe in good agreement with the 541+89 
PNe observed by iReid & Parker! (|2006). Since most of these PNe come from single or non- 
interacting binary stars in our model, this means that most such stars produce PNe contrary to 
the "Binary Hypothesis" which suggests that binary interaction is required to produce a PN. 

Key words: binaries: close - planetary nebulae: general - stars: late-type. 



1 INTRODUCTION 

In the LMC and SMC, red giant variables follow at least five 
Pe riod-Luminosity re lations, originally named sequences A to 
E dWood etal.lll999h . Sequence A, B and C are radial pulsa- 
tion sequences and sequence D, which is parallel to sequences 
A, B and C but of longer period, is stil l of unknown physi- 



A, a and I but ot l onger pe riod, is stil l ot unknown physi- 
cal origin dWood. Olivier. & Kawaled |2004 iNicholls et al.ll2009l 



iNie. Zhang. & Jiand2010h . Sequence E stars follow a loose Period- 
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Luminosity relation between sequences C and D, but they extend 
to lower luminosities. Sequence E stars are widel y accepted to be 
binaries jWood et al.lll999l;ISoszvnski et al.ll2004l) . The orbital pe- 
riods of the sequence E red giants are typically 100-600 days and 
the full light curve amplitudes are less than 0.6 mag in the MA- 
CHO red (M R ) band. They lie on the first giant branch in the case 
of low-mass stars, or the equivale nt pre-core helium burning giant 
phase of intermediate mass stars l lWood et"aL 1999 | ). Some of the 
sequence E star s also extend to the AGB l Soszvnski et alj|2004t 
Soszyriski 2007). As a whole, sequence E stars make up approxi- 
mately 0.5-2% of the RGB and AGB stars in the Large Magellanic 
Cloud (LMC). 
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Sequence E vari ables show ellipsoidal light variations 
dSoszvhski etafl |2004> . In these systems, the primary star has 
evolved to the red giant phase and substantially filled its Roche 
lobe. Because the orbital separation of the system is similar to the 
size of the red giant, the red giant is distorted from spherical sym- 
metry to an ellipsoid or pear-like shape by tidal interactions with 
its unseen companion. It is the orbital rotation of this shape that 
gives rise to the ellipsoidal light variations. About 7% of sequence 
E va riables show eclipses in addition to the ellipsoidal light varia- 
tions dSoszvhski etal]|2004) . 

The first radial velocity curve s of se quence E stars were pre- 
sented in iNicholls. Wood, & Cionil j201 d) and from their observa- 
tions, the mean full velocity amplitude is about 43 km s -1 . This is 
consistent with the amplitudes expected for close red giant binaries 
with roughly solar-mass components. The velocity variation of the 
ellipsoidal variable is regular, dominated by its orbital motion. For 
every single cycle of the velocity variation there are two cycles of 
light variation, the latter being caused by the change in apparent 
surface area. 

In past and ongoing radial velocity studies of sequence E stars 
in the LMC, we have obtained spectra for 110 systems. In none of 
these systems do we see the emission lines characteristic of symbi- 
otic stars. This means that the companion is unlikely to be a white 
dwarf or neutron star. It is presumably a main-sequence star so that 
the red giant is usually the primary star in the binary system. 

As the Roche lobe of the sequence E red variables is already 
substantially filled, further expansion of the red giant as it evolves 
is expected to soon cause Roche lobe overflow. The binary system 
will suffer a common envelope (CE) event (except in the rare case 
that the red giant has previously lost sufficient mass that mass trans- 
fer is stable). A CE event will lead to the ejection of the red giant 
envelope and the termination of RGB or AGB evolution. If the two 
stellar cores do not merge in the CE event, a luminous compact 
close binary system will be left. In the case of a CE event occur- 
ring on the RGB, the compact system will be a low mass white 
dwarf binary; while on the AGB, the remnant stellar core will be the 
central star of a planetary nebula (PN). If the two stellar cores do 
merge the y should form a rapidly rotating red giant of the FK Co- 
mae class dWebbinkll976l ; lBopp & Stencelll98ll ; lBopp & Rucinskil 
Il98ll). The fraction o f the rapidly rotating red giants is about 2% 
dCarlberg et ai1l2oT 1), although many of these stars may be ti dally 
locked binaries rather than single stars. 

The most common morphology of PNe is nonspherical and 
there are several theories attempting to explain the underlying 
mechanisms responsible for the shaping. Among all the mech- 
anisms, the binary hypothesis is most favoured (Bond & Livid 



199d:lYungelson. Tutukov. & Livid 1993k ISokerlll997l : lBondl200d : 
Ziilstral2007f De Marcdl2009l) and it has been confirmed by obser- 



vations in some cases. In fact, the shape of PNe with close binary 
central stars is definitely observed to be asymmetric. It has been 
suggested that binary int eraction plays a major role in the forma - 
tion and shaping o f many (IBond & Livid 199d:llben& Livid 19931) . 
or even most PNe dBondl200d ; lDe Marcdl2009l) . 

At the present time, the actual fraction of PNe with close bi- 
nary central stars is not clear. The short periods of the known cen- 
tral stars of binary PNe mean that these stars have been through a 
CE phase since the companion currently lies well within the radius 
of the precursor star when it was a red giant. It seems fairly clear 
that the brighter sequence E stars are indeed the immediate precur- 
sors of the known binary PNe central stars. At lower luminosities, 
the low mass white dwarf binaries, which are p roduced from RGB 
binaries (e.g. |Paczvnski|[l976l : IWebbinkll 19841) . should be descen- 



dants of sequence E stars. It is our aim to use the observed fraction 
of sequence E stars among all red giants to estimate the fraction 
of PNe with close binary central stars. We also use the sequence 
E fraction and the overall population of RGB stars in the LMC to 
estimate the birthrate and population of all PNe there, along with 
the birthrate of low mass white dwarf binaries. 

In a study usin g MACHO l i ght cu rves of stars in a central 
part of the LMC bar, IWood et alj j 1999b found that -0.5% of the 
red giants in the top one magnitude of the RGB were detected as 
sequen ce E variables in the M R band. Similarly, ISoszvnski et al.1 
d2004l) and lSoszvnskH d2007l) using OGLE II data found that 1-2% 
of red giants were detected as sequence E variables. These numbers 
provide the calibrating input for the fractions of close binary PNe 
and low mass helium white dwarfs produced by CE events. Simul- 
taneously, with this calibration, we can also estimate the fraction of 
red giants that are in wider binaries which reach the AGB tip with- 
out Roche lobe filling, and the fraction of PNe that are descended 
from single stars. 



2 THE SIMULATION MODEL 

The prediction of the evolutionary fate of binary and single red gi- 
ants is made by a Monte Carlo simulation. One million red giant 
binary systems are generated using orbital element distributions de- 
rived from observations of stars in the solar vicinity. The evolution- 
ary fate of these binary systems is then examined. Single stars are 
added to the total stellar population in order to reproduce the ob- 
served fraction of sequence E stars on the red giant branch. Details 
of the standard simulation model are described in Section |2~T1 Be- 
sides the standard model, we adjust the inputs to see how dependent 
the results are on the model assumptions. Results obtained by vary- 
ing one or multiple parameters are described in Section|4] Since we 
will be utilizing the observed sequence E population in the LMC, 
our model uses inputs from LMC sources when available. 

2.1 The standard model 

2.1.1 Generating the binary system 

(1) The initial mass of the primary star 
Initially we require the red giant to be the primary star in the binary 
systems. The primary mass is drawn from the initial mass distribu- 
tion /j(m, t) defined by 



dN = h(m, t)dmdt 



(1) 



where him, t) is the number of stars born per unit mass per unit time. 
To obtain h(m, t), we consider the Initial Mass Function (IMF) as 
well as the star formation history. Let g(m) be the IMF (assumed 
independent of time) and let fit) be the total star formation rate at 
time t, then h(m, t) is expressed as: 

h(m,t) = g(m)xf(t) . (2) 

Here, g(m) is assumed to follow the ISalpeterl d 19551) power law: 
g(m) = m~ 2,35 . An observa tionally-based appro ximation to the 
LMC star formation history dBertelli et al.lll992h has a star burst 
which starts ~4 Gyrs ago and ceases ~0.5 Gyrs ago, with the ratio 
of burst to quiescent star formation rate 2?burst being 10. Thus, /(f) 
is expressed as: 

( , t < 0.5 

fit) = \ fltar,, x b , 0.5 < t < 4 , 

[ b, t>4 
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where i?burst = 10, b is a constant, and t is in G yrs from the presen t 
time. According to the evolutionary tracks of iGirardi et al. 
appropriate for the LMC (Z=0.008), after converting the ages of 
red giants to initial masses, we have: 



Mm)) = k(m) 




m > 3.0 
3.0 > m > 1.3 
m < 1.3 



(3) 



where m is in solar masses. 

With the initial mass distribution, the cumulative probability dis- 
tribution of the mass, H(m), is expressed by: 



I g(m) X k(m)dm 

Jmo 



(4) 



where mo corresponds to the oldest stars currently on the giant 
branch. 

As we only consider red giant stars, the mass range should be that 
of RG B and AGB stars. The evolutionary tracks of IGirardi et al.l 
( feOOOh show that the minimum mass of low-mass stars in the LMC 
that are currently on the first giant branch is about 0.9 M Q . Stars 
with m < 0.9 M Q have not yet evolved to the RGB since such stars 
have a main sequence lifetime longer than the time since the first 
stars formed. Hence, in Equation l(4} we set m to 0.9 M Q . The up- 
per limiting mass for RGB stars which develop electron degenerate 
helium cores on the first ascent of giant branch, is about 1.85 M Q . 
Stars with mass larger than this do not reach the high luminosi- 
ties of the observed sequence E stars in the LMC until they enter 
the early-AGB (EAGB). As the LMC star formation is assumed to 
cease 0.5 Gyrs ago, the maximum mass for AGB stars is 3.0 M Q . In 
the Monte Carlo simulation, we adjust 'b' so that H(3.0) = 1, and 
draw the primary mass from the normalized H(m). 

(2) The initial mass ratio 

We define q = m 2 lm l as the mass ratio, where m\ is the mass of the 
primary star and mi is the mass of the secondary star. 

The initial mass ratio is drawn from the distribution given in 
iDuquennov & Ma yor (1991) who give orbital element distributions 
for binary stars in the solar vicinity. The mass ratio distribution fol- 
lows a Gaussian-type relation, with a peak at 0.23. The reason we 
use the distribution from the solar vicinity (as well as the orbital pe- 
riod distribution described below) is that the mass ratio distribution 
is unknown observationally in the LMC. 

Given that we now know both masses in the binary system, we 
can decide which component is 'currently' being seen as the red 
giant. This is done by randomly assigning the primary or the sec- 
ondary as the current red giant with probabilities that are propor- 
tional to the star formation rate at the time the system would have 
been born if the primary or secondary, respectively, was the current 
red giant. In the case where the secondary is the current red giant, 
the full evolution of the primary is followed before the evolution of 
the secondary is considered. 

(3) The orbital period 

T he initial orbital period (P) is drawn from the distribution given 
in IDuquennov & MavoJ dl99ll) . The period distribution follows a 
Gaussian-type relation, with a peak of log.P=4.8, in unit of days. 

Note that some bin aries with long periods are too wide to influ- 
ence the shape of PNe. lSoker, Rappaport. & Harpazl l l 1998b suggest 
P < 500-2000 years if a wide binary is to influence a PN shape. 
iMastrodemos & Morrisl ( 1 19991) similarly give a gravitational focus- 
ing fraction ar/ oc greater than 0.03 for a binary to produce a shell 
more distorted than "quasi-spherical". If we apply this to low-mass 
stars with masses of ~1.5 M Q and wind velocities of ~10 km s" 1 , 
we find P < several hundred years is required. Based on the above, 



we adopt 500 years as the upper limiting period P max for a wide 
binary to influence the shape of a PN. Hereinafter, we refer to bi- 
naries with P > P max as 'wide binaries' and binaries with P < P max 
but with periods large enough to avoid Roche lobe overflow on the 
RGB and AGB as 'intermediate period' binaries. 

(4) Eccentricity and stellar separation 

About 90% of the observed light curves of sequenc e E stars indicate 
zero or small eccentricity dSoszvnski et alj|2004l) . so we assume 
zero eccentricity. Given the primary mass, mass ratio and orbital 
period, and assumed zero eccentricity, we can calculate the separa- 
tion of the stellar centres a via the equation of orbital motion. 

(5) Orbital inclination 

The orbital inclination ( is obtained assuming a random orientation 
of the orbital pole so that dN = sin idi. 



2.1.2 Luminosity limits for RGB and AGB stars 

According to observa t ional studies dFrogel. Cohen. & Perssonl 
ll983l : IWood et al.ll 19991: iKiss & Bedding||2003l) and the the evolu- 
tionary tracks of IGirardi et al.ld2000t) for a typical LMC metallicity 
of Z=0.008, the bolometric magnitude for the RGB tip (TRGB) is 
close to -3.6 mag. Thus, we set Mt, i(TRGB) = -3.6 as the maxi- 
mum luminosity of sequence E stars on the RGB. In our simulation 
we compare the observed and simulated ratio of sequence E stars 
to all red giants at luminosities corresponding to the brightest one 
magnitude of the RGB, i.e. -2.6 > M bol > -3.6. We consider the 
evolution of stars once they become brighter than the base of the 
RGB. 

AGB stars extend in luminosity up to the point at which a 
superwind rapidly terminates the AGB evolution. The luminosity 
at the tip of the AGB, M bo i(TAGB), is calculated from: 



Mboi(TAGB) = -5.1 - 0.6(m,. - 1.55) 



(5) 



Here, the zero point -5.1 is obtained from the maximum luminos- 
ity for AGB stars in the r ich intermediate-age LMC cluster NGC 
1978 dKamath et"ai]|2010h . In t his cluster, the mass of the O-rich 
stars on the RGB is -1.55 M Q dKamath et"aHl20ld) . In Equation 
(|5j, m-ifj is the initial mass of the red giant. The slope 0.6 is the 
AGB tip mag nitude change per solar mass for stars in the range of 
1.0-3 .0 M Q dFrogel. Mould. & Blanco] 1 1990l : IVassiliadis & Wood 
1 1994b . Binary stars that evolve to the AGB tip without filling their 
Roche lobes are assumed to eject their envelopes rapidly by a su- 
perwind, in the same way as a single star. 



2.1.3 Mass loss from the red giants 

We consider stellar mas s loss for all ev olving stars. We use the 
empirical formulation by iReimersI d 1 9751) to calculate the mass loss 
rate, but with the rate multiplied by a parameter n which is set equal 
to 0.33 dlben & Renzindll983l:lLebzelter & Woodll2005h . 

According to Tout & Eggletonl 1 19881) , the mass loss rate may 
be tidally enhanced by a factor of: 

R 1 

1 +Bxmin[( — ) 6 , — ] , 
R L ' 2 6 ' 



where R is the radius of the primary sta r , R L i s the equivalent ra- 
dius of the Roche lobe. iTout & Eggletonl dl988l) suggest B= 10000. 
However, it is unlikely that such enhanced mass loss rates are re- 
alistic for a wide range of binary systems (see the references in 
Section[4]l, so we treat 'B' as a variable parameter. In our standard 
model it is set to zero. 

In order to calculate the amount of mass lost per magnitude of 
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evolution up the giant branch, and to calculate the lifetimes of the 
sequence E stars, we need the evolution rate dM bol /dt. For low-mass 
stars wit h m < 1.85 Ma on th e RGB, the evolutionary rate in the 
tracks of iBertelli etail j2008b is well approximated (in mag/Myr) 
by 



dM bL 
dt 



-7.58 x irr 



■ 2.96 x lfr 4 — + 4.73 x 1CT 8 ( — | 



On the EAGB, where the helium burning shell alone supplies most 
of the star's luminosity, the evolution rate for all initial masses con- 
sidered here (0.9-3.0 M Q ) is approximated by 



dM hL 
dt 



1.15 + 0.75M bl 



(7) 



On the AGB {M bo \ < -3.6), the average evolu tion rate of the full 
AGB calculations o flVassiliadis & Wood! d 19941) for LMC metallic- 
ity is approximated for all masses in the range 0.9-3.0 M Q by 



dM b , 
dt 



1/(0.456 -0.085M,; ol ) 



(8) 



Although these formulae are derived from evolutionary models for 
single stars, they also apply for red giants in binary systems since 
the nuclear evolution of the red giant c ore is independent of the 
conditions in the conve ctive envelope dRefsdal & Weigertlll970l : 
I Wood &Faulkneill 19861) . 

Many binary stars in the simulation reach luminosities on the 
AGB where a superwind causes AGB termination, just as happens 
for a single star. This termination is assumed to occur over a brief 
interval and it is treated by specifying an AGB termination lumi- 
nosity (see Sec tion l2, 1,21 ). 

Since the primary star loses mass in a wind, the secondary star 
can accrete some of t he material as it orbits t hrough it. We include 
wind accretion by the Bondi & Hovle ll 19441) mechanism, as given 
by equation (6) of iHurlev. Tout. & Pols! d2002l) with aw = 3/2 and 
pw = 1 /2. The orbital evolution resulting from mass loss and ac- 
cretion is calculated using equation (20) in IHurlev et al.l d2002l) . 



2.1.4 Tidal effects 

When a red giant substantially fills its Roche lobe, tides induced 
in the red giant by the orbiting companion cause the conversion of 
orbital angular momentum into spin angular momentum of the red 
giant. We treat this process using the equations in Section 2.3.1 of 
IHurlev et al.l d2002h . In particular, the rate of change of spin angu- 
lar mo mentum is given by equations (26) and (35) of IHurlev et al.l 
d2002h . with the spin angular momentum change being extracted 
from orbital angular momentum of the binary system. We also in- 
clude the loss of spin angular momentum from the red giant b y 
stellar wind mass loss using equation (1 1) o f lHurlev et al1d2002h . 



2.1.5 Light variation of an ellipsoidal variable 

For an ellipsoidal binary system containing a red giant, the ampli- 
tude of the light variation depends on the fractional filling of the 
Roche lobe by the red giant, the mass ratio and the orbital inclina- 
tion. 

Rather than usi ng an analytic approximation to the ellipsoidal 
light variations (e.g. [Morris I ll985l) , we use the program nightfaliQ 
to model the observed ellipsoidal light variations of partial-Roche 



1 http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall.html 



lobe filling systems as a function of binary parameters. In these 
models, the effective temperature of the red giant is set to 4000 K, 
as most of our sequence E stars are K or M type stars. A blackbody 
flux distribution is assumed. Since the companions of the sequence 
E stars are expected to be mostly main sequence stars, the radius 
and luminosity of the secondary star are set to small values so they 
do not affect the light variation . As noted in Section 12.1.11 we as- 
sume zero eccentricity. We then use nightfall to create light curves 
for ellipsoidal variables with a range of values for the Roche lobe 
radius filling factor /, mass ratio q and inclination ;'. Then full light 
amplitudes in the M R and / bands are derived, and a fit to the input 
parameters is made. We find the full light variation amplitude in 
M R band is well approximated by: 



AM R = (0.22 1/ 4 + 0.005) x (1.44956? 025 



0.44956) x sirti , (9) 



for 0.5 < / < 0.9, 0.1 < q < 1.5, < i < § . We also find the 
relation between AM S and AI is well approximated by: 



Al = 0.87 x AM R 



(10) 



where AI is the full light curve amplitude in the / band. The lines 
in Fig.Q]show fits to the factors in equations (|9) and ( 110) . 



2.1.6 Radii during the ellipsoidal phase of the red giant 
evolution 

Given a minimum amplitude for detectable light variations (see 
Section |2.1.9t , the minimum radius filling factor f miu of a red gi- 
ant in a binary system with detectable ellipsoidal light variation 
can be determined from Equations (|9) and dlOt after q and i have 
been selected. The corresponding minimum red giant radius R miu is 
given by 



^min — X f m 



(ID 



where R L is the equivalent radius of the Roche lobe dEggletor] 
Il983l) : 



0.49^ 



-2/3 



Rl = 

a 0.6q- 2 ' 3 + ln(l + q- [ ' 3 ) ' 

The maximum stellar radius R mia , corresponding to the Roche lobe 
being filled, will be 



R n 



R, 



(12) 



2.1.7 Roche lobe overflow 

The further expansion of a Roche lobe-filling red giant as it evolves 
will cause Roche lobe overflow. Then, provided the mass ratio q = 
m 1 lm l is smaller than the critical mass ratio q a \i = (m 2 /m i ) C y it for 
the unstable mass transfer (see below), the binary components will 
come into contact and undergo a CE event. 

In binaries that undergo a CE event, there are two possible out- 
comes - a close binary system or a single coalesced star. Follow- 
ing w idely used procedures (e.g. see Section 2.7.1 of IHurlev et al .1 
|2002|) . we assume that the orbital energy dissipated in the interac- 
tion process goes into overcoming the binding energy of the red 
giant envelope which is thus lost. If neither the remnant red giant 
core nor the companion main sequence star fills its Roche lobe at 
the end of this process, a close binary containing a white dwarf (the 
red giant core) and a main sequence star will be left. However, if the 
main sequence star fills its Roche lobe, then Roche lobe overflow 
will cause the main sequence star to mix with the common enve- 
lope. The red giant core thus gains H-rich material and the result 
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1.44956q" ""'-0.44956 




Figure 1. Fits of AMr to the factors in Equation l|9) (top three panels) and 
the relation between Al and AMr (bottom panel). In the first and third pan- 
els, blue diamonds denote the values of AMr when q = 0.5, green squares 
denote the values of AMr when q = 1.0, and the red triangles denote the 
values of AMr when q = 1.5. Line fits according to Equation l|9) are also 
shown, blue dotted lines correspond to blue diamonds, green long dashed 
lines correspond to green squares, and the red dot-long dashed lines corre- 
spond to red triangles. In other panels, solid circles denote the values ob- 
tained from the nightfall program and the corresponding black lines are the 
fits. 



will be a single, merged red giant whose further evolution will pro- 
duce a PN with a single-star nucleus. Note that as the red giant core 
is much more compact than the main sequence star, it is unlikely 
that this core will fill its Roche lobe. 

The above p ossibilities are treate d in our model using equa- 
tions (69)-(77) of iHurlev et all (120021) . w ith A = . 5 and a cE = 1 
in our standard model . Accor ding to Ide Kool dl990L Il992l): 
lYungelson et alj J 1993b and Irian. Podsiadlowski. & Eggletonl 
dl995d) . a CE = 1 gives the best fit to observations (we also try 
lower values - see Section[4). We assume the mass of the red giant 
core is given by the formulae L/L Q = 59250 x (M C /M G - 0.48) 
for AGB stars and L/L = 5.3 x 10 5 (M c /M o ) 6 - 667 for RGB stars, 
which are the fits to the core mass-luminosity relati on on the red 

[ iant branch shown by the evolutionary tracks of iBertelli et al.1 
200g|) at metallicity Z=0.008. 

To check when stable m ass transfer ra t her th an a CE event oc- 
curs, we use equation (57) of IHurlev et al.l ( 120021) to calculate q ail . 
Note that this equation assumes conservative mass transfer which 
we also assume. Except for small envelope masses (less than 20% 
of the red giant mass), the critical mass ratio (tni/mi)^ given by 



the above equation is greater than 1.0 and very few of red giants 
end up with stable mass transfer in our standard model. Very high 
mass loss rates or different values of q ait can give larger numbers 
of red giants with stable mass transfer (see Section|4]where we run 
models with a lower value for q ali ). 

In binaries that undergo stable mass transfer, the red giant 
evolves closely filling its Roche lobe and transferring matter to 
its companion. Roche lobe filling is maintained in the models by 
enforcing a large mass transfer rate when the red giant radius ex- 
ceed its Roche lob e radius. This is done using equation (2) of 
IChen & HanU2008h . 

2.1.8 Effective temperature and bolometric magnitude 

Low-mass stars evolve through the RGB and later ascend the 
EAGB and AGB. When stars enter the EAGB and AGB they are 
oxygen rich. However, after several thermal pulses, which lead to 
dredge-up of carbon to the stellar atmosphere, the O-rich stars are 
converted to C-rich stars. This leads to a decrease in effective tem- 
perature r cff . Owing to this, T c fs values are calculated separately for 
these two sub- types. 

We use the intermediate age LMC globular cluster NGC 1978 
as a template for LMC red giant properties. According to the data 
of iKamath et~al] d201Cl) . the giant branch slope in the HR dia- 
gram for O-rich stars is given by d l og T^/dM^ = 0.03 5. More- 
over, from the evolutionary tracks of iGirardi et al.l (feOOOl) . we find 
dlogTeff/dm = 0.03 (with m in solar masses). The effective tem- 
perature of O-rich stars is thus given by 



log r eff = log r cff . + 0.035M bol + 0.03m, 



(13) 



To find the zero point log T c g$ for the effective temperature, we 
fit Equation J 1 3 b to the observed HR diag ram (Fig. [2^ for LMC 
sequence E stars in ISoszvnski et al.l ( 12004 ). The effective temper- 
atures of the observed sequence E stars are calc ulated by convert- 
ing the ir (J - K) to r eff using the tr ansforms inlHoudashelt et al.l 
d2000dlbl) . We adopt E(B - V)=0.08 dKeller & WoodkOOfj). E{T- 
K) = 0.52£(fi - V) and A(K) = 0.35E(B - V) dRieke & Lebofskvl 
Il985l) . The bolometric correction BC K is calculat ed from (J - K) 
with the transforms in iHoudashelt et al ] J2000allbh . The J and K 
magn itudes were obtained from the 2MASS catalogue dCutri et al .1 
120031) . Fitting the RGB given by Equation (13) for masses between 
0.9 M Q and 1.85 M s to the densest part of the observed RGB in 
Fig.[2]gives logr cff , = 3.648. Note that 0.9 and 1.85 M Q are the 
lower and upper mass limits, respectively, for low mass RGB stars 
(see Section l2.1.1t . 

There are substantial numbers of sequence E red giants in Fig. 
|2]which are brighter than M bol = -3 and which have log T cff about 
0.08 larger than that of the low mass RGB stars. These red giants 
are intermediate mass stars (M > 1.85 M G ) with non-degenerate 
He cores evolving to the point of He core burning or perhaps on 
the EAGB. Other sequence E red giants with large deviations in 
log T c g from the main locus of the RGB in Fig. [21 such as those in 
the lower right of the figure, presumably have photometry contam- 
inated by spatially coincident field stars of the LMC or the Galaxy 
or they may be high mass loss rate AGB stars with thick obscuring 
circumstellar shells. 

Because C-rich stars are often surrounded by substantial dust 
shells that cause an unknown amount of reddening, it is not possi- 
ble to use observed photometric colors to derive the photospheric 
temperature. We therefore use red giant models to obtain the po- 
sition of the C-rich star giant branch. For C-rich stars, the slope 
of the giant branch in the HR diagram is obtained from Table 5 
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Figure 2. The HR diagram of the sequence E stars. The red points are 
the observation data, while the black lines are the red giant branch fits for 
masses of 0.9 M G and 1.85 M Q . Typical observational error bars are shown 
at three magnitude levels. The maximum luminosity of each black line cor- 
responds to the transition luminosity from O-rich to C-rich stars. The lumi- 
nosity of the RGB tip is also marked. 




Light full amplitude (mag) 

Figure 3. The relative distribution of full light amplitude for ellipsoidal 
variables on the top one magnitude of the RGB. Grey lines are estimates 
from the Monte Carlo si mulation and the black lines are the observed 
OGLE II distribution from ISoszvnski et alj J2004 I) (botto m panel) and the 
observed MACHO distribution from lWood et alj d 19991) (top panel). The 
histograms in the top panel apply for the MACHO red band while the his- 
tograms in the bottom panel apply for the OGLE II / band. 



in iKamath et al] feOlCh . The models in this table were computed 
to fit the O-rich giant branch in NGC 1978, as well as allowing 
for C/O ratios to increase from unity at luminosities brighter than 
the transition luminosity for O to C-rich stars. The models predict 
that dlog r cff /rfM bo i =0.125 for C-rich stars of ma ss ~1.55 M,. If 
we ad opt the distant modulus to the LMC of 18.54 dKeller & Wood 
I2006I) . the transition luminosity from the O-rich to C-rich stars is 

M bol = -3.57 - 0.6m! , 

based on the data of iFrogel et alj dl990l) and IVassiliadis & Woo3 
dl994l) . Matching the effective temperature for O-rich and C-rich 
stars at the transition luminosity, the effective temperature of C-rich 
stars becomes 



log T ei 



3.969 + 0.084m, + 0.125M b , 



(14) 



With L = 4no-R 2 T* s and log L/L Q = -0.4(M bol - 4.75), we 
can derive the bolometric magnitude for O-rich and C-rich stars 
as a function of stellar radius R from Equations d!3t and <l 14b . For 
O-rich stars, 



M bo , = 4.362 - 0.222m, - 3.7041ogfl 

and for C-rich stars, 

M bol = 1.189 - 0.373m, - 2.2221ogfl 



(15) 



(16) 



If we substitute R with R m [ n or R mdx , then we get the minimum 
or maximum bolometric magnitude of the red giants with just- 
detectable light variations or with Roche lobes that are just full, 
respectively. We note that the formulae given above for T eff are 
derived from single star observations. However, these formulae 
should also give good results for red giants in binary systems as 
stellar radii are not greatly affected by the presence of a binary 
companion l |Paczvhskj||l97l[ and references therein). 



2.1.9 The MACHO and OGLE comparison data 

As noted in the introduction, IWood et aL d 19991) and 
ISoszvnski et al .1 d2004l) found that the fraction of red giants 
which have detectable ellipsoidal variabili ty is ~0.5% and 1-2% , 
respectively. The higher fraction found by ISoszvnski et alj d2004l) 
is due to the smaller photometric errors of the OGLE observations 
and hence the ability of OGLE to detect smaller ellipsoidal 
variations than MACHO. The fraction of stars that show ellipsoidal 
variability is an important input to our calculations, so we recom- 
pute these fractions accurately. We also investigate the amplitude 
required for detectability of ellipsoidal variability in each of these 
two studies. 

Using the MACHO data of I Wood et all d 19991) and selecting 
the top one magnitude of the RGB, we find that the fraction of red 
giants with detectable ellipsoidal variability is 0.61%. Similarly, us- 
ing the OGLE II ellipsoidal var i ables on the top one magnitude of 
the RGB from lSoszvnski et"al] 1 2004[) and the c orresponding total 
red giant population from llJdalski et alj d2000l) . we find an ellip- 
soidal fraction of 1.33%. 

The distributions of full light amplitude for the MACHO and 
OGLE II ellipsoidal variables are shown in Fig. [3] It is clear that 
the OGLE II observations detect variability at considerably lower 
amplitude than the MACHO observations. Also shown in Fig.[3]are 
the corresponding distributions of ellipsoidal light amplitudes on 
the top one magnitude of the RGB computed from our Monte Carlo 
model. We take the detectability limits for ellipsoidal variability 
to be the light amplitude at which the observed number of stars 
has fallen to 50% of the estimated true number. These limits are 
0.055 magnitudes for the M R band data and 0.019 magnitudes for 
the OGLE / band data. We use these detectability limits and the 
observed relative fractions of ellipsoidal variables to calibrate our 
Monte Carlo model. 
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2.1.10 A summary of evolutionary outcomes 

We aim to predict the fractions of binary and single red giant stars 
that suffer each of the various possible evolutionary fates. The pos- 
sible evolutionary outcomes are itemized below. 

(1) The star reaches the AGB tip without filling its Roche lobe. 
Single stars will follow this path and produce a single PN. Wide and 
intermediate period binaries will also follow this path, producing a 
PN with its nucleus in a wide or intermediate period binary. 

(2) The star fills its Roche lobe on the AGB but above the RGB 
tip. There are three possible outcomes. 

(a) A CE event leads to the formation of a close binary con- 
sisting of an AGB star core and a secondary star. In this case, 
evolution of the AGB star core leads to creation of a PN with a 
close binary nucleus. 

(b) A CE event leads to a merger of the two stars with the 
secondary star merging into the red giant envelope. The merged 
star will then evolve to the AGB tip as a single star. 

(c) Stable mass transfer occurs leading to the formation of 
an intermediate period binary. If all the mass lost from the red 
giant was transferred to the secondary, no PN would be visible. 
However, there is likely to be some mass lost from the system 
and we assume a PN will be seen in the rare cases that stable 
mass transfer occurs. 

(3) The star fills its Roche lobe below the RGB tip. For stars 
with initial masses m < 1.85 M Q , the Roche lobe filling will occur 
on the RGB. However, for stars with m > 1.85 M G , this will occur 
on the EAGB. There are three possible outcomes of such CE events. 

(a) A CE event leads to the formation of a close binary con- 
sisting of a red giant core and a secondary star. In this case, evo- 
lution of the red giant core to high effective temperatures is too 
slow to l ead to creation of a PN (see the post-RGB evolution 
tracks of Drie be et al. I GUI) but note the comments at the be- 
ginning of Section[3j the cores of EAGB stars also evolve slowly 
as they are burning He in a shell). We will generally call these 
stars post-RGB stars even though some of them come from the 
EAGB. Theoretical considerations suggest that the outcome of a 
CE event in a close binary on the RGB (or EAGB) is likely to 
be a post-CE binary consisting of a low-mass He-core (or CO- 
core) white dwarf and a basically unaltered low-mass secondary 
JPaczvnslri|[l976l : IWebbmld[l984l) . We suggest that some Pop- 
ulation II Cepheids, UU Herculis variables and RV Tauri stars 
with luminosities below the RGB tip are also formed in this way. 
Most, if not all, RV Tauri stars are binary system s that have mid- 
infrared dust emission and circumbinary disks dde Ruvter et al.l 
120051 ; I Lloyd Evans & van Winckell2007l) , as expected from ter- 
minating red giant evolution on the RGB or AGB by Roche lobe 
filling. Observationally, some of the Population II Cepheids, UU 
Herculis variables and RV Tauri sta rs in the LMC have luminosi- 
ties below the RGB tip luminosity jAlcock et al J 19981) . 

(b) A CE event leads to a merger of the two stars with the 
secondary star merging into the red giant envelope. The merged 
star will then evolve to the AGB tip as a single star. 

(c) Stable mass transfer occurs leading to the formation of 
an intermediate period binary. Some of the binary Population II 
Cepheids, UU Herculis variables and RV Tauri stars with dust 
disks could be produced in this way. 



2.1.11 Simulation normalization 

We use the well determined fractions of the sequence E stars on 
the top one mag nitude of the RGB (Section |2~L9l Wood's 0.61% 
and Soszynski's 1.33%) to normalize our calculations. Given our 
adopted IMF, period distribution and mass ratio distributions, some 
binary stars will have parameters such that the Roche lobe fills or 
nearly fills at these luminosities (-2.6 > M bol > -3.6), and the star 
will exhibit detectable ellipsoidal light variations characteristic of 
sequence E stars over a luminosity sub-interval of one magnitude. 
The lifetime as a sequence E star, corresponding to this luminosity 
sub-interval, is calculated by using the evolutionary rate described 
in Section [2. 1.31 Some of the other binaries will have their evo- 
lution terminated by interaction before reaching M bol = -2.6 and 
some will evolve through the full interval without showing any de- 
tectable light variations. In the latter case, the lifetime is also cal- 
culated using the evolutionary rate described in Section |2. 1.31 We 
calculate and sum the sequence E lifetimes of all RGB and EAGB 
stars in the interval -2.6 > M bol > -3.6 and we also calculate and 
sum the total lifetimes of all red giants as they pass through this 
interval. The ratio of the summed sequence E lifetimes of those bi- 
naries showing detectable ellipsoidal variations to the total summed 
lifetimes of all binaries gives the apparent fraction of sequence E 
stars. This apparent fraction due to binaries alone is always greater 
than the observed fraction. Therefore, after we have considered 10 6 
binary stars in the Monte Carlo simulation, we add single stars to 
the population of stars, assuming they evolve through the full in- 
terval -2.6 > M bo | > -3.6, until the predicted and observed se- 
quence E fractions match. Note that we are implicitly assuming 
here that the ratio of binary to single stars in the LMC is a free 
parameter and that it can differ fro m the ratio estimated f o r stars 
in the solar vicinity. According to iDuquennov & Mayor! d 199 lh 
the percentage of sol ar-type stars that are single is ~33% while 
iRaghavan et al .1 (120 ldh, who derive a perio d distribution quite sim- 
ilar to that of lDuquennov & Mayor! ( 1 199 lh . estimate the percentage 
of single solar-type stars to be ~56%. 

We emphasize that our main purpose in this study is to use 
the observed sequence E star population, which has initial orbital 
periods of ~ 100-600 days, to estimate the population of red giants 
whose evolution is affected by binary interactions. These interact- 
ing stars typically have initial orbital periods of <3000 days. The 
observed period distribution of binary stars (IDuquennov & MavoJ 
1 1 99 it IRaghavan et al .120101) is very broad with a peak at ~ 10 5 days 
and the stars we are primarily concerned with occupy only a small 
part of the overall period distribution. Our predictions related to 
red giant evolution should thus be relatively insensitive to the over- 
all period distribution (provided it is smooth). However, estimates 
of the numbers of stars with periods >>3000 days, and the number 
of single stars, will depend more sensitively on the overall period 
distribution. 



3 RESULTS FOR THE STANDARD MODEL 

In this section, we present the results for our standard model using 
both the Wood and Soszynski frequencies for sequence E stars. The 
models produce the relative birthrates of stars that terminate their 
AGB or RGB evolution. The numbers in Table[TJgive these relative 
birthrates. 

Throughout this paper we make some assumptions. Firstly, 
we assume that all stars that leave the AGB produce a PN. Sec- 
ondly, when we compare ratios of populations of different types 
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of PNe (for example close binaries PNe and all PNe) we assume 
that the PN lifetimes is independent of the AGB termination pro- 
cess (CE event or wind ) although this may not be the case (e.g. 
iMoe & De Marcel 20061) . Similar considerations apply for the post- 
RGB stars. It should be kept in mind that we are really comparing 
birth rates. 

An example of how evolution rates could be different for post- 
CE stars and single or non-interacting binaries is provided by the 
class of post- AGB (perhaps post-RGB) star that consists of a binary 
system with a dusty circumbinary disk, an or bital period of a few 
hundred days or more, and a large eccentricity (Ivan Winckell2007h . 
These stars could result from AGB or RGB termination by stable 
mass transfer. Alternatively, their evolutionary path could be caused 
by the high eccentricity of their orbits. The disks around these stars 
may provide a reservoir of hydrogen that can be accreted back onto 
the high luminosity star to fuel nuclear burning for an extended pe- 
riod, at the same time keeping the hydrogen envelope thick enough 
to maintain a relatively low T cff which is insufficient to excite a PN. 

At the other extreme is the possibility that a common enve- 
lope event in a RGB star can remove enough of the hydrogen above 
the burning shell that the star comes out of the CE event with r eff 
large enough to excite a PN almost immediately dMoe & De Marco] 
12006b . In our modelling, we assume that this does not happen and 
that the evolution of post-RGB stars is so slow that the circumstel- 
lar shell has dispersed before r c ff reaches values high enough to 
produce a PN (see Section l2. 1.101 item 3a). 

3.1 Fractions of binary planetary nebula nuclei and 
post-RGB stars 

With the standard inputs and the evolutionary scenarios described 
above, we examine which binary systems form close binary PNe, 
intermediate period binary PNe, wide binary PNe, single binary 
PNe, and post-RGB binaries. The predicted PNe fractions of the 
standard model are listed in the first two lines of TableQ] In the total 
population of planetary nebula nuclei (PNNe), close binaries make 
up 9% or 7%, intermediate period binaries make up 27% or 23%, 
wide binaries make up 55% or 46%, and single stars make up 3% or 
19%, using Wood's (0.61%) or Soszynski's (1.33%) frequency, re- 
spectively. The ratio of post-RGB star births to close binary PNNe 
births is ~50%. Of those binaries undergoing a CE event on the 
AGB, we find ~1.4% suffer a merger while the remainder produce 
a close binary. Due to the relatively higher binding energy of stars 
on the RGB, the frequency of mergers happening on the RGB is 

much higher than that on the AGB. Below M bo \ 1.5, essentially 

all red giants with a main sequence companion merge when Roche 
lobe overflow occurs. However, if the companion is a white dwarf, 
a double degenerate binary can be produced down to low lumi- 
nosities. Our standard model predicts that merged stars, which are 
likely to be rapidly rotating red giants, would make up about 5% 
of the red giants on the RGB above M bol 1.5. L ower on the gi- 
ant bra nch, the merged fraction will be less than 5%. ICarlberg et al] 
d201 ll) estimate the observed fraction of rapidly rotating red giants 
to be -2%. 

Our results can be compared to the binary evolution models 
of lHanetalHl995al) . Their models 4 and 1 1 have input parameters 
similar to our standard model except that our initial period distribu- 
tion is different from theirs (they assume a constant number of bi- 
naries per interval of log a, where a is the binary separation). In our 
standard model, we find th at about 8% of PNNe are close binaries 
whereas lHan et al] dl995a]) predict that about 4-5% of PNNe are 
close binaries (see the column "CE Ejection (AGB)" in Table 1 of 



lHan et all l Tl995al) and column "Binary PNNe / CE" in this paper). 
The reason for our higher fraction can be traced to the different ini- 
tial period distributions. Binary systems that undergo CE ejections 
on the AGB have initial orbital periods of about 230-1400 days or 
separations of about 230-760 R Q (assuming the binary component 
masses are ~1.5 M Q ). All binary systems with initial periods (or 
separations) greater than these values will produce PNe through 
wind mass loss, as will single stars. The number of binaries with 
initial periods (or separations) that lead to a CE event on the AGB, 
relative to the number of all binaries plus single stars that produce 
PNe, is rou ghly twice as large in our standard model as in models 
4 and 1 1 of lHan et al] dl995al ). This explains the higher fraction of 
close binary PNNe produced in our standard model. (We note that 
our models 15 and 16, described in Section HI have s imilar period 
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distributions to models 4 and 1 1 o f lHan et alTTl995 j) . Both sets of 
models predict that about 4—5% of PNNe are close binaries indicat- 
ing similar outcomes for similar input physics.) 

Finally, we note that the fractions of close binary PNNe given 
by lHan et all Jl995id) are somewhat arbitrary since these fraction 
depends on the maximum binary separation allowed for their bi- 
nary systems (they adopt, without giving justification, a maximum 
separation of 5.75 x 10 6 R G ). They also make the assumption that 
there are no single stars. In our models, the fraction of stars that are 
single or in wide or intermediate period binaries is observationally 
constrained by the fraction of stars showing ellipsoidal variability 
on the RGB. 

We now examine observational estimates of the fraction of 
close binary PNe. Searches carried out by Bond and collabora- 
tors for close binary central stars of PNe using photometric vari- 
ability te chniques obtained a fraction of 10-15% for close bi- 
nary PNe dBond & GraueJl987MBond. Ciardullo & Meakesll992t 
lBondlll994l.l2000h . However, the survey bias is not well understood 
since the PNe samples were monitored over 30 years with many dif- 
ferent observing campaigns dDe Marco. Hillwig. & Smithl2008 ). A 
new survey for close binary central stars by [Miszalski et al. I d2009l) 
got a fraction of 12-21% for close binary PNe. This survey is 
more efficient and was carried out in a relatively uniform manner, 
by searching for periodic photometric variability of homogeneous 
PNe samples. However, their 12-21% is not a definitive fraction for 
close binary central stars. For example, the four close binary central 
stars PHR 1744-3 355, PHR 1801-2718 PHR 1804-2645 and PHR 
1804-2913 in the iMiszalski et al] d2009l) s ample of 22 have been 
subsequently questioned and excluded by [Miszalski et al] d201ll) 
which would reduce their close binary fraction to 10-16%. On the 
other hand, accounting for the effect of orbital inclination on the 
light curve amplitude, the derived close binary fraction could in- 
crease. 

Our standard model prediction of ~7-9% for close binary 
PNe, is lower than the fractions estimated from the observations 
mentioned above. However, our simulation is strongly constrained 
by the observed fractions of sequence E stars. In addition, our 
model reproduces the observed light amplitude, period and velocity 
amplitude distributions of sequence E stars as well as the fraction of 
low mass He white dwarfs in the total white dwarf population (see 
Section [3, 3. 21 . This suggests that our estimated fraction of PNe 
with close binary central stars is reliable. 



Table 1. Model parameters and relative birthrates of red giant progeny 
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No a C E 


^max 


^burst 


IMF 


Pdist 


delist 


B 


q 






Single PNNe 




Binary PNNe 


Binaries 




Single 


Wide 


MgAGB 


MgRGB 


IntP 


Stable 


CE 


CE 


Stable 


PRGB/PP 


Standard model 




































1 1.0 


500 


10 


S1955 


DM1991 


DM1991 





H57 


w 


3.31 


54.90 


0.12 


5.71 


27.20 


0.01 


8.74 


100.00 


0.00 


0.0448 


















s 


18.52 


46.24 


0.11 


4.83 


22.97 


0.01 


7.31 


100.00 


0.00 


0.0375 


One parameter variation models 
































L U.O 


son 


J u 


*s 1 QSS 

C3 1 7JJ 


DM1991 


DM1991 





H57 


w 


1.67 


55.00 


0.55 


7.28 


27.21 


0.01 


8.28 


100.00 


0.00 


0.0275 


















s 


17.18 


46.28 


0.46 


6.10 


22.95 


0.01 


7.01 


100.00 


0.00 


0.0233 


J yJ.J 


son 


i n 
J u 




DM1991 


DM1991 





H57 


w 


0.00 


54.99 


1.83 


8.81 


27.31 


0.01 


7.05 


100.00 


0.00 


0.0101 


















s 


16.62 


45.97 


1.50 


7.29 


22.73 


0.01 


5.88 


100.00 


0.00 


0.0085 






i n 
J u 


Q 1 Q S S 


DM1991 


DM1991 





H57 


w 


19.77 


0.00 


0.13 


5.89 


65.38 


0.01 


8.81 


100.00 


0.00 


0.0461 


















s 


33.00 


0.00 


0.11 


4.86 


54.59 


0.01 


7.44 


100.00 


0.00 


0.0388 


5 1 .0 


1 00 


1 


S1955 


DM1991 


DM1991 





H57 


w 


0.00 


68.71 


0.13 


5.80 


16.62 


0.01 


8.72 


100.00 


0.00 


0.0441 


















s 


14.55 


59.29 


0.11 


4.80 


13.93 


0.01 


7.31 


100.00 


0.00 


0.0370 


a in 

O L.V 




i 
i 


c 1 Q s s 


DM1991 


DM1991 





H57 


w 


3.08 


55.46 


0.07 


5.71 


27.99 


0.08 


7.61 


100.00 


0.00 


0.0650 


















s 


18.67 


46.68 


0.06 


4.78 


23.35 


0.07 


6.41 


100.00 


0.00 


0.0546 


7 1.0 


500 


10 


Reid02 


DM1991 


DM1991 





H57 


w 


1.10 


56.35 


0.19 


5.59 


27.63 


0.01 


9.12 


100.00 


0.00 


0.0407 


















s 


15.48 


48.04 


0.17 


4.75 


23.72 


0.01 


7.83 


100.00 


0.00 


0.0342 


8 1.0 


500 


10 


SI 955 


Flat 


DM1991 





H57 


w 


62.15 


3.91 


0.11 


10.45 


16.37 


0.01 


6.99 


100.00 


0.00 


0.0480 


















s 


66.87 


3.42 


0.11 


9.10 


14.39 


0.01 


6.11 


100.00 


0.00 


0.0419 


9 1.0 


500 


10 


SI 955 


DM1991 


n(q) = q- QA 





H57 


w 


11.40 


41.73 


0.20 


6.61 


30.57 


0.02 


9.47 


100.00 


0.00 


0.0533 


















s 


23.55 


36.11 


0.17 


5.68 


26.32 


0.02 


8.15 


100.00 


0.00 


0.0461 


10 1.0 


500 


10 


SI 955 


DM1991 


n(q) = 1 





H57 


w 


27.55 


34.43 


0.09 


4.93 


25.37 


0.03 


7.61 


100.00 


0.00 


0.0447 


















s 


35.73 


30.60 


0.08 


4.37 


22.50 


0.02 


6.70 


100.00 


0.00 


0.0398 


11 1.0 


500 


10 


SI 955 


DM1991 


n(q) = q 





H57 


w 


37.00 


30.58 


0.03 


3.69 


22.62 


0.04 


6.05 


100.00 


0.00 


0.0439 


















s 


41.04 


28.62 


0.02 


3.43 


21.13 


0.04 


5.71 


100.00 


0.00 


0.0411 


12 1.0 


500 


10 


SI 955 


DM1991 


DM1991 


500 


H57 


w 


8.16 


52.63 


0.10 


5.33 


27.32 


0.73 


5.74 


99.30 


0.70 


0.0478 


















s 


19.65 


45.97 


0.09 


4.71 


23.96 


0.63 


4.98 


99.28 


0.72 


0.0420 


13 1.0 


500 


10 


SI 955 


DM1991 


DM1991 





0.8 
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23.07 


43.35 


0.10 


4.42 


21.79 


0.72 


6.55 


92.18 


7.82 


0.0367 


















s 


26.68 


41.35 


0.10 


4.26 


20.74 


0.68 


6.19 


92.02 


7.98 


0.0351 


Multiple parameter variation models 
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R2010 
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16 38 


42.63 


0.10 
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28.20 


0.03 
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0.00 


0450 


















s 


26.15 


37.63 


0.08 


4.23 


24.92 


0.02 


6.96 


100.00 


0.00 


0.0399 


15 1.0 


100 


1 


MS 1979 


Flat 


»(<?) = 1 





H57 


w 


63.43 


8.96 


0.05 


11.00 


10.79 


0.10 


5.67 


100.00 


0.00 


0.0720 


















s 


66.63 


8.13 


0.05 


10.12 


9.73 


0.12 


5.23 


100.00 


0.00 


0.0650 


16 1.0 


100 


1 


MS 1979 


Flat 


n(q) = q 





H57 


w 


71.36 


7.44 


0.01 


7.96 


8.91 


0.13 


4.19 


100.00 


0.00 


0.0654 


















s 


71.99 


7.29 


0.02 


7.71 


8.66 


0.13 


4.21 


100.00 


0.00 


0.0630 


17 1.0 


500 


10 


SI 955 


DM1991 


DM1991 


500 


0.8 


w 


26.15 


42.00 


0.08 


4.24 


22.16 


1.59 


3.77 


86.64 


13.36 


0.0420 


















s 


28.44 


40.62 


0.08 


4.17 


21.44 


1.56 


3.68 


86.06 


13.94 


0.0396 



Model parameters: acE is the orbital energy transfer efficiency (see Section [2. 1.71 ; P max is the ma ximum initial p eriod (in years ) for which a bin ary is assumed t o influence the shape of PNe; Rburst is the star 
formation rate burst factor (se e Secti on |2Tjj; I MF is the assumed initial mass function ('S195 5' fromfsalpeter 19551 'Reid02' from lReid et al.l2002l 'MS 1979' from lMiller & Scalolll979h : Pdist is the initial period 
distribution ('DM1991' from lDuquennov & Mayor 199lJ 'R2010' from lRaghavan et al.|[2oToL 'Flat' means dN oc dlogP); q^t is the initial mass ratio distribution n(q) cc dN/dq; B i s the parameter in t he tidally 
enhanced mass loss rate formula (see Section. [2.1. 3) ; q C nt = fe/fllcrii is the critical mass ratio above which Roche lobe overflow is stable ('H57' is the value given by equation (57) oflHurlev et al. 2 0021) . 



Model results: Lines containing 'w'/'s' use Wood's/Soszyriski's sequence E frequency to normalize the results. The numbers in columns 11-19 give the relative birth rates for PNNe and post-RGB stars (as a 
percentage). The sum of all PNNe birthrates is 100%, as is the sum of the post-RGB star birthrates. The ratio of the total post-RGB star birthrate to the total of the PNNe birth rate is given column 20. PNNe 
types are: stars that are born single (Single); wide binaries which have P > P max and which do no influence the shape of their PNe (Wide); binary stars that merge on the AGB (MgAGB) or the RGB (MgRGB); 
intermediate period binaries that never fill their Roche lobes but which have P < P mllx (IntP); binary stars that fill their Roche lobes and then undergo stable mass transfer (Stable); and binary stars that fill their 
Roche lobes and then undergo a CE event leaving a close binary (CE). Post-RGB are: those that result from CE events (CE); those that result from stable mass transfer (Stable). 
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100 




Figure 4. The orbital period distribution of sequence E stars. The red 
solid l ine denotes the distribution from the observations oflSosz yriski et al] 
120041) . the green dotted line denotes the distribution from the model using 
Wood's frequency, and the blue short dashed line denotes the distribution 
from the model using Soszyhski's frequency. 



3.2 Properties of sequence E stars 

To further test our simulation, we predict the binary properties of 
sequence E stars and compare them with observations where pos- 
sible. The main available data are the photometric variations from 
MACHO and OGL E experiments, which provide a probe into the 
period distribution. iNicholls et al.l 1 120 101) provide radial velocity 
observations of sequence E stars, although only 1 1 binary systems 
were studied. 



3.2.1 Period distribution 

In Fig. [4] we present the orbital period distribution of sequence E 
stars, predicted with Wood's and Soszyhski's frequencies. The or- 
bital period contribution of a single binary to the distribution is the 
average of the periods at onset and end of detectable sequence E 
variability on the RGB, weighted by the sequence E lifetime for 
the system: note that the orbital period typically changes by less 
than 3% during the sequence E phase (due to mass loss and tidal 
interaction). We als o show the distribution o f the observed periods, 
using the sample o f ISoszvhski etal] ( 120041) . We only consider the 
sequence E stars on the top one magnitude of the RGB, so the stars 
in the observed period distribution are selected based on their lumi- 
nosity. Fig. [4] shows that the simulation reproduces the observation 
very well: not only are the total numbers of sequence E stars re- 
produced, as required, but their period distribution is reproduced as 
well. We also note that the shape of the light amplitude distribution 
is reproduced well by the model (Fig. [3}. 

3.2.2 Full velocity amplitude distribution 

In Fig. [5] we show the full velocity amplitude dist ributions from 
our m odel and the radial velocity observations by INicholls et al.l 
feoioh . The expected full velocity amplitude concentrates in the 
range of ~5 to 50 km s -1 , with a peak at ~22.5 km s -1 . It shows a 
reasonable agreement with the observations, although the observed 
number is small. 
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Figure 5. The full velocity amplitude distribution of sequen ce E stars. The 
red so lid line denotes the distribution from the observation of Nicholls et al. 
l201fj|) with 1 1 stars, the green dotted line denotes the distribution from the 
model using Wood's frequency, and the blue short dashed line denotes the 
distribution from the model using Soszyhski's frequency. The distribution 
peaks are normalized to 1 . 




Figure 6. The mass distribution of the sequence E stars. The green dotted 
line denotes the distribution from the model when using Wood's frequency 
while the blue short dashed line denotes the distribution from the model 
when using Soszyhski's frequency. 



3.2.3 Mass distribution 

In Fig. [6] we show the mass distribution of the sequence E stars. 
The mass range is from 0.85-3.0 M Q , with the LMC star burst be- 
ginning at 1.3 M Q . The lower limiting mass of 0.85 M Q is set by 
the mass of the oldest stars and stellar wind mass loss. There are a 
large number of stars in the range of 1.3-1.85 M G , where the lat- 
ter mass is the maximum mass for stars that ascend the RGB with 
an electron degenerate core. The more massive stars in Fig. [6] are 
EAGB stars at luminosities corresponding to the top magnitude of 
the RGB. 
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Figure 7. The same as Fig.|6]but for the mass ratio. 
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Figure 8. The same as Fig.|6]but for the orbital separation. 



falls to zero at 2 or 4 Myr. Longer lifetimes are predicted in the 
Soszynski case than the Wood case. The reason is that OGLE II has 
a more sensitive detection limit for photometric variations (~0.02 
mag) than MACHO (~0.05 mag) (see Fig.[3}, so Soszynski's stars 
are detected initially with less-filled Roche lobes and it then takes 
longer before the lobe is filled and RGB evolution is terminated. 
The average lifetime for sequence E stars is ~0.95 Myr, when the 
ellip soidal variation amplitude i s more than 0.02 magnitudes. 

lMikoiaiewsk3 12007. 2011) claims that in symbiotic stars, the 
red giant radius is often close to half the Roche lobe radius yet 
ellipsoidal light variations are seen in these stars with an ampli- 
tude that suggests the Roche lobe is almost filled. If this claim is 
correct, then it suggests some unidentified source of light variation 
that simulates an ellipsoidal variation. The sequence E stars are not 
as extreme as the symbiotic stars where the red giant is usually a 
semi-regular or Mira variable with a large mass loss rate and the 
companion is a hot accreting white dwarf whose radiation may af- 
fect the facing surface of the red giant. We therefore do not expect 
our ellipsoidal light curve calculations for sequence E stars to be af- 
fected by an unidentified source of ellipsoidal-like light variation. 
Nevertheless we investigate how such an extra source would alter 
our results if it did exist. To simulate the extra light amplitude sug- 
gested in symbiotic stars, we replaced the Roche lobe filling factor 
/ in Equation l|9} by 2f. The effect of this is to make detectable 
ellipsoidal variability occur when the radius of the red giant is half 
the usual value. Note that this does not affect the evolution of any 
star, it simply means that ellipsoidal variability is detectable in red 
giants with wider orbits than in our usual calculations. This in turn 
causes a larger fraction of binaries to appear as sequence E stars. 
Now, in order to reproduce the observed fractions of sequence E 
stars, a much larger population of single stars is required than in 
our usual models. For example, in our standard model, about 10% 
of stars are single and 90% are in binaries whereas in the models 
with the modified Equation (|9) about 80% of stars are single and 
20% are in binaries. This means that the fraction of all stars that un- 
dergo CE events or mergers is reduced by a factor of 90/20 = 4.5. 
In particular, the fraction of close binary PNNe is reduced from 
~8% to ~ 1 .8%. (We also note that in our modified models the ellip- 
soidal light amplitudes of stars nearly filling their Roche lobes are 
increased by a factor of 16 giving rise to ellipsoidal light curve am- 
plitudes up to 3 magnitudes, which are never observed (see Fig.f3j. 
Of course, the modified Equation ([9} has no physical basis so we 
cannot really say what the ellipsoidal amplitudes should be.) 



3.2.4 Mass ratio distribution 

In Fig. [Vj we show the mass ratio distribution of sequence E stars. 
The mass ratio is less than 1.1, with a peak at ~0.5. Mass loss 
from the red giant is the reason that some systems have mass ra- 
tio slightly greater than 1 . 

3.2.5 Orbital separation distribution 

In Fig. [8] we show the orbital separation distribution of sequence E 
stars. The separation is in the range of ~ 100-600 Rq, with a narrow 
peak at -200 R Q . 

3.2.6 Ellipsoidal variability lifetime 

We show the predicted ellipsoidal variability lifetime of sequence 
E stars in Fig. [9] The distribution starts at a peak of ~0.2 Myr due 
to the contribution of EAGB stars with m > 1.85 M Q , and then 



3.3 Properties of binary PNNe and post-RGB stars 

3.3.1 Luminosity distribution 

In Fig. [10] we show the luminosity distributions of the outcomes of 
the binary red giants according to the standard model. As expected, 
intermediate period binary PNNe have luminosities of the AGB tip, 
close binary PNNe have AGB luminosities above the RGB tip, and 
EAGB and RGB binaries undergoing a CE event but not merging 
have luminosities below the RGB tip. Mergers occur preferentially 
at luminosities lower than the luminosities of CE events that lead 
to binaries. The RGB binaries formed at M bo] > -1.5 are double 
degenerates. In these systems, the small radius of the companion 
white dwarf allows them to avoid a merger. 
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Figure 9. The same as Fig.|6]but for the lifetime of the ellipsoidal variabil- 
ity. 
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Figure 10. Luminosity distributions of binary red giants at the termination 
of their red giant evolution. Panel (a): intermediate period binary PNNe 
with P < 500 years. Panel (b): close binary PNNe (blue) and merged 
AGB stars (green). Panel (c): EAGB binaries (red) and merged EAGB 
stars (cyan/grey). Panel (d): RGB binaries (red) and merged RGB stars 
(cyan/grey). In all panels, dotted lines denote the distributions from the 
model when using Wood's frequency, short dashed lines denote the dis- 
tributions from the model when using Soszyhski's frequency. 



3.3.2 Mass distribution 

In Fig.[TT] we show the mass distributions of the binary red giants 
at their termination luminosities. As expected, stars in intermedi- 
ate period binaries lose a significant amount of their mass via the 
Reimers-like stellar wind. Mergers occur preferentially at lower lu- 
minosity so such systems have the least mass loss. 

At the end of the superwind phase or a CE event, a red giant 
has lost its whole envelope, leaving a He (or CO) core as a rem- 
nant. In Fig. [12] we present the mass distribution after the red giant 
envelope has been lost. The mass of the close binary PNNe is in the 
range of ~0.50-0.80 M Q , with most of the stars concentrated be- 
tween 0.50-0.65 M Q . For the post-RGB binaries, the mass is lower, 
with a range of ~0. 2-0.45 M . The post-EAGB CO-cores concen- 
trate around 0.50 M G . The mass distributions of intermediate period 
binary PNNe are also shown in Fig. [12] The masses are in the range 
of ~0.55-0.85 M Q , with a strong peak at ~0.6 M . 

Mass determination of 86 PNNe by Stasiriska & Tvlendal 

dl990h shows a mass distribution similar to our prediction: the mass 
range is 0.55-0.75 M Q , an d most of the PNNe con centrate between 
0.55-0.65 M . Similarly, IZhang & Kwokl d 19931) determined the 
mass of 303 PNNe, reporting a mass distribution from 0.55-0.85 
Mo, with a narrow peak at 0.6 M H . Th e population synth esis model 
results o f lYungelson et al.l d 19931) and lHan et all dl995al) give mass 
distributions of close binary PNNe (with q-ce = 1) that are similar 
to ours. 

For the post-RGB stars, the predicted masses are con- 
sistent with observ ations of low-mass white d warf bina- 
ries. For example. iLiebert. Bergeron. & Holbergl (12005b and 
iRebassa-Mansergas et alj i 201ll) determined the mass distribu- 
tion of local white dwarfs and found a distribution with a peak 
at 0.4 Mq as well as a more dominant peak around 0.6 M Q 
which corresponds to CO white dwarfs produced by evolution 
off the AGB. The low mass peak is produced predominantly by 
CE events on the RGB, although some of t hese stars (<30%) 
may be single dKilic. Brown. & McLeodl20ld : iBrown et al]|201ll ; 
IRebassa-Mansergas et alj|201ll . and references therein) which re- 
quires an exotic type of evolution such as an unusually high wind 
mass loss on the RGB, or ejection of the envelope by a giant planet, 
or mergers of t wo He white dwarfs ( we d o not find any such events). 
The st udies of ILiebert et al.l d2005 ) and IRebassa-Mansergas et all 
d201 ll) suggest that the birth rate of low mass He white dwarfs is 
about 4-10% of the total white dwarf birth rate. Our standard model 
predicts this relative birth rate to be 9.6%, in good agreement with 
the observations. Our predicted mass distributio n is also consis- 
tent w ith th e observati o ns and the simulations of [Yungelson et all 
dl993l) and lHan et alj dl995j) . Note that the He white dwarfs 
with masses from 0.2-0.3 M G in Fig. [12] are nearly all double 
degenerates and they form about 19% of the He white dwarfs. 
Some He white dwarfs are known to be in double degenerate sys- 
tems dMarsh. Dhillon. & Duckl 19951) although the observed frac - 



tion is very poorly known dKilic. Stanek. & Pinsonneaulj 120071) . 
Our models predict that only ~0. 1 % of sequence E stars will have a 
white dwarf companion (on the way to producing a double degener- 
ate), consistent with our finding mentioned in the introduction that 
there is no evidence for degenerate companions in a sample of 1 10 
sequence E stars. 

Finally, we note that our adopted star formation history means 
that there are no stars with initial masses greater than 3.0M o . Even 
if we were to allow recent star formation by (say) using f(t) = b for 
t < 0.5, the relative number of stars with M > 3.0 M Q would be tiny 
because of the combined effects of the IMF and the star formation 
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Figure 11. Mass distributions of the binary red giants according to their 
evolutionary fates. Thin lines are for the initial mass distributions. Thick 
lines are for the mass distribution at the tennination luminosity: the top 
panel shows the mass for intermediate period binaries at the AGB tip and 
just before the beginning of the superwind phase; the other panels show the 
mass at the start of the CE event. Colors are as in Fig. 1101 



history. We would not simply be extending the plots in Figs. II llto 
higher masses according to the IMF. We would be extending the 
plots at a level reduced by a factor of 10 in height since the stars 
with mass 1.3-3.0 M G were produced during the burst. 



3.3.3 Period distribution 

In Fig. Q~3] we present the orbital period distribution of the close 
binary PNNe and post-RGB binaries. The orbital periods of the 
close binary PNNe have a range from 0.01-1000 days, and a peak 
at ~4 days. The orbital periods of post-RGB binaries are between 
0.01-10 days, and have a narrow peak at ~0.3 days. 

Orbital period distrib ution s of close binary P NNe predicted 
bv lYungelsoneUu1 ll993l) and lHan et alj dl995ah have a similar 
distribution to ours, showing a period range of few hours to 1000 
days. Observ ationally, periods of close binary PNNe are all less 
than 16 days teondl2000l : iDe Marcdl200^ ; lMiszalski et alj|200gh 
with the current period distributions peaking at or below P ~ 1 day. 
This peak is at shorter periods than predicted by the model (the 
blue curve in Fig.ll3t. The absence of the long period PNNe in the 
current sample could be due to the observational b ias against long 
period stars: in variability searches such as that of iMiszalski et al.1 
d2009h . short periods in close binaries are needed to get the de- 
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Figure 12. Mass distributions of intermediate period binary PNNe (black), 
close binary PNNe (blue), post-RGB and post-EAGB binaries (cyan/grey) 
and double degenerate secondaries (red). Note that the cyan curves include 
the double degenerate component. 
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Figure 13. Orbital period distributions of the objects shown in Fig. [12] Line 
types are the same as in that figure. 



tectable light amplitudes from ellipsoidal and reflection effects; 
and in radial velocity searches, close binaries with short periods 
are needed to get the high velocity variations needed for these hot 
compact stars which have relatively broad lines. Alternatively, the 
current observed period distribution and the model period distribu- 
tions can be brought into better agreement by using a lower value 
of acE (see Section|4](. 

The orbital period distribution of local binary white d warf s 
(e.g. iRebassa-Mansergas et al.1 120081 ; IZorotovic et al.1 1201 Oh also 
peaks at periods less than or near 1 day. The model prediction given 
by the cyan line in Fig. [T3] has a peak around 1 day but more stars 
with periods longer than 1 day than the current observational esti- 
mates. 
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Figure 14. Mass ratio distributions of the objects shown in Fig. [12] Line 
types are the same as in that figure. 



of iReid &~P arker 2006) so we adopt 7.5 arcsec as a typical maxi- 
mum diameter for LMC PNe. If we use a distance modulus to the 
LMC of 18.54 feller & Woocj|2006h and a PN expan sion velocity 
v exp o f 34.5 km s _1 (the median for the 94 LMC PNe in lDopita et all 
1 19881) then we find tpn = 2.6 x 10 4 years. Using the above num- 
be rs, we predict that the number of PNe in the region surveyed 
by IReid & Parked j2006h will be mrgb+eagb X TmlT\ma& = 548 
PNe. The estimated numb er of PNe is in good a greement with the 
541±89 PNe observed bv lReid & Parker! l2006ri . Since more than 
90% of these PNe come from single or non-interacting binary stars 
in our model, this means that most such stars produce a PN. This 
is contrary to th e "Binary Hypothesis' ' dMoe & De Marco| [2006: 
|Pe Marcoll2009l) which suggests that binary interaction is required 
to produce a PN. 

There are a small number of AGB star terminations that may 
not produce a detectable PN. Most of the single or non-interacting 
binary stars in our model terminate their AGB evolution by a super- 
wind. However, some have their evolution terminated by the low 
mass loss rate Reimers wind. These are the stars with the lowest 
initial mass 0.9 < M/M s < 0.92 which make up only 0.7% of 
the AGB ter minations. It is these stars that may not produce a de- 
tectable PN l lSoker & SubagboOSh . 



3.3.4 Mass ratio distribution 



In Fig.[l4]we show the mass ratio distributions of close PNNe and 
ost-RG B binaries. This is similar to the prediction of lHan et al.l 
1995ah . 



3.4 The population of PNe in the LMC 

We use our model to predict the number of PNe i n the inner 25 
square degrees of the LMC searched for PNe by IReid & Parkerl 
d2006h . As part of our modelling procedure to estimate the apparent 
fraction of sequence E stars in the top one magnitude of the RGB 
(-2.6 > Mb i > -3.6), we compute the average lifetime Ti mag in 
this magnitude interval for all RGB and EAGB stars which evolve 
into this interval. We find ri mag = 2.77 x 10 s years for the standard 
model, with a deviation from this value of less than 3% for other 
models. The average lifetime is determined by the evolution rate 
and variation in the average lifetime is caused by the small fraction 
of stars whose lives are terminated by CE e vents in this interval. 
From the SAGE catalog jMeixner et al ■I2006I) . we obtained the near 
and mid-inf rared photomet r ic obs ervations of all stars in the area 
searched bv lReid & Parkerl d2006l) . Using the position of the giant 
branch in the (J,J— [3.6]) diagram, we selected all the stars (RGB 
and EAGB) in a parallelogram corresponding to the top one mag- 
nitude of the RGB. The parallelogram has sides J= 13. 25, J=14.25, 
J-[3.6 ]=3.828 - 0.17J and J-[3.6]=3.328 - 0.17J (see lBlum et all 
for SAGE colour-magnitude diagrams for the whole LMC). 
We found the total number of stars hrgb+eagb = 58356 and es- 
timate an error of less than 5% in this number. The uncertainty 
is due to scattering of stars into and out of the parallelogram in 
the (J,J-[3.6]) diagram because of observational errors and con- 
tamination by stars at the red end of core helium burning loops. 
Confusion is not a problem for these relativ ely bright stars (see 
iNikolaev & Weinberdl2000tlBlum et alj2006h . 

The final quantity we need in order to estimate the number of 
PNe in the inner 25 square degrees of the LMC is the average life- 
time t pn of LMC PNe. We note that there is a steep drop in the num- 
ber of LMC PNe with diameters between 7 and 8 arcsec (see Fig. 16 



4 RESULTS FOR MODELS WITH VARIED INPUT 
ASSUMPTIONS 

In order to investigate how dependent the results are on the model 
parameters (e.g. ckce, IMF, period and mass ratio distribution), as 
well as to compare the results with other population synthesis cal- 
culations, we have run models with adjusted input parameters. 
Models have been made varying one, two or three of the param- 
eters described above and the results are shown in Table [T] We dis- 
cuss the effects on each model of varying the model parameters. 
For selected models, distributions of final masses and periods are 
plotted in Figs. 1 151 and 1161 which are similar to Figs. 1121 and 1131 
respectively. 

Models 2 and 3: the energy efficiency parameter «ce- 
The value of the energy efficiency parameter acE is controver- 
sial and there is no consensus for its value. As noted in Sec- 
tion l2.1.71 a value of 1.0 seems consistent with some observational 
constraints. However, smaller values of acE have also been derived 
along with estimates of dependency on the parame t ers of the binary 
dPolitano & Weilerl l2007r. IZqrotovic et al.l l201ol : be Marco etaD 
1201 it iDavis. Kolb. & Kniggdl201ll) . To investigate the effect of 
changes to a CE on our modelling results, the standard value of 1.0 
was changed to 0.6 and 0.3 in models 2 and 3, respectively. As ex- 
pected, decreasing q-ce causes more mergers to occur at the expense 
of fewer CE events (Table [T). Reducing a C E from 1 to 0.3 changes 
the fraction of close binary PN produced from ~ 8% to ~ 6.5%. The 
results for model 3 are shown in Figs. [l5]and[T6] The most inter- 
esting thing about these plots is that the periods of the close binary 
PNNe now peak just shortward of 1 day, in much better agreement 
with the current observed period distribution. A lower value of ctce 
is a way of removing the discrepancy between the observed period 
distribution and the period distribution predicted by our standard 
model. 

Models 4 and 5: The upper limiting period P max . 

P m has been changed from the standard value of 500 ye a rs to in - 

finity or 100 yrs. The latter value was used by |Hanetal1 l ll995ah . 
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Figure 15. Mass distributions of intermediate period binary PNNe (black), 
close binary PNNe (blue), post-RGB and post-EAGB binaries (cyan/grey) 
and double degenerate secondaries (red) in models 1 (thick solid lines), 3 
(thick dashed lines), 1 1 (thin solid lines) and 12 (thin dashed lines). Results 
are shown only for the models using Wood's sequence E frequency. 
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Figure 16. Orbital period distributions of the objects shown in Fig.[l5] Line 
types are the same as in that figure. 



The main effect of this parameter is to shift stars between the inter- 
mediate period, wide and single categories. 

Model 6: The star burst ratio /f hurst' 

In this model, the star burst ratio /Jhurst m Section |2.1.1| was changed 
to 1 .0 to si mulate a constant star f ormation rate, which is often used 
in models jde Koo<ll99d Il992l ; lYungelson et"ai]| 19931 ; lHan et al.l 
Il995ah IPolitano & Weileill2007t) . The change produces more stars 
of mass less than 1.3 M relative to higher mass stars. The effect 
is a small reduction in the number of CE events on the AGB and 
an increase in the number of CE events on the RGB. This is prob- 
ably due to the smaller initial orbital separations of the lower mass 



binaries, given that the same period distribution is assumed at all 
masses. 

Model 7 -.The IMF. 

Here t he Salpeter IMF was changed to the lReid. Gizis. & Hawlevl 
d2002h power law (with an exponent of -1.3). The effect on the 
birthrate ratios is small, with slightly more CE events on the AGB 
and fewer on the RGB. 

Model 8: The period distribution P ji st . 

Here a flat period distribution with dN cc dlogP was used. This 
period distribution has been commonly used (see Ide Kooll 1 1990l . 
1 19921 : iHanetaill 1995 j) . The period distribution has a large influ- 
ence on the relative birthrates of single, wide and intermediate pe- 
riod PNNe. Since the flat period distribution has far fewer stars 
born at long periods relative to the periods of the sequence E stars, 
it is necessary to add many more single stars than in the standard 
model in order the get the fraction of sequence E stars on the top 
magnitude of the RGB to agree with the observational value. The 
fraction of stars undergoing CE events is slightly lower than in the 
standard case. This model is a good indicator of the fact that using 
sequence E stars to normalize our results keeps the fraction of CE 
events in our simulations fairly constant when model parameters 
are changed. 

Models 9 to 11: Mass ratio distributions q^ st . 
Initial mass ratio distributions n(q) = Ccf with a = -0.4, 
or 1 were tried. The first exponent was that found by 
iKouwenhoven et all J2007I ) for main-sequence A and B stars (M ~ 
2-15 M Q ) in the Sco OB2 association. The higher the fraction of 
high mass companions, the fewer the number of mergers, as might 
be expected from the higher orbital energy available to eject the red 
giant envelope. The higher the fraction of high mass companions 
also gives a smaller percentage of CE events in our model. This is 
probably because, for a given fraction of sequence E stars, a higher 
mass companion can give a detectable ellipsoidal light variation at 
a larger orbital separation, reducing the likelihood of a CE event 
occurring. The final mass and period distributions for model 1 1 are 
shown in Figs. \T5\ and [16] These plots clearly show an increase 
in the fraction of double degenerates produced when the number of 
high mass companions is increased. The periods of the close binary 
PNNe are also increased slightly. 



Model 12: Tidally enhanced mass loss. 

The tidally enhanced mass loss rate o f tTout&EggletonlJ 19881) was 
tried, with the free para meter 'fi' set to 500, although its actual 



value is very uncertain ( Han et al 


Il995bl; 


Hani 19981: Soker et al. 


1 19981; iKarakas, Tout, & Lattanzic 


l200d; 


Frankowski & Tvlenda 


1200 ll; iHurlev et alj 20021). The main effect is an increase in the 



number of stars undergoing stable mass transfer (because of the in- 
crease in m 2 /mi resulting from mass loss from the red giant) and a 
decrease in the number of CE events. The final mass and period dis- 
tributions for model 12 are shown in Figs.[l5]and[T6] The reduced 
number of CE events on the AGB is clearly evident in the figures. 
There is also an extension of the period distribution of close binary 
PNNe to higher values, further increasing the discrepancy between 
the observed and predicted period distributions. The enhanced mass 
loss also significantly increased the minimum period of intermedi- 
ate period PNNe. 

Models 13 and 17: Critical mass ratio q alt . 

The critical mass ratio q crit = (m 2 lm{) ait , above which mass trans- 
fer from a Roche lobe filling red giant is stable, has some uncer- 



16 J.D.Nieetal. 



taintv. lChen&HarJ (2008) investigated the value of q ait appropri- 
ate for mass transfer from red giants under various conditions. For 
the more luminous stars and for mass transfer efficiencies more than 
0.5, the value of q clit is mostly similar to that used here in our stan- 
da rd model. However, values as low as (m2/mi) cl j t = 0.8 are found 
bv lChen&Hanlte008h in some cases. In order to see how depen- 
dent our models are on q crit , we ran a simulation with q cnt = 0.8. We 
also ran a simulation with q aA = 0.8 and '6' = 500 to see how the 
extra reduction in the red giant mass from the enhanced mass loss 
would affect the number of stars undergoing stable mass transfer. 
As expected, the number of stars undergoing stable mass transfer is 
greatly increased, especially on the RGB. The fraction of binaries 
undergoing CE events is decreased because of the alternative evolu- 
tionary path for Roche lobe filling stars. Changing q alt to 0.8 makes 
no significant difference to the mass and period distributions. 



Model 14: The \Raeliavan et al\ \201&) distributions. 
Here we use the p eriod and mass ratio distributions of 
iRaghavan et which ca n be considered as update d 

versions of the dist ributions of iDuquennov & Mavorl dl99lh . 
Raghavanetai](l2010h lave a period distribution similar to that of 



Duquennov & Mavoi J 1 99 lh but they have a flat mass ratio distri- 



bution. The results are very similar to those of model 10 which also 
has a flat mass ratio distribution. 



Models 15 and 16: The\Miller &■ Scald \l97<X) IMF. 
Here we use the lMiller & Scald dl979h lognormal IMF appropriate 
for low-mass stars, combined with a flat period distribution, and 
n( q) = 1 or q. Our model 15 has similar parameters to model 4 
of lHan et alj Jl995iJ) . As with model 8, the flat period distribution 
leads to a high single star fraction in our simulation. The combined 
effect of all parameters is to decrease the number of CE events on 
the AGB and increase the number on the RGB. 



We estimate that about one third of PNe contain binaries ca- 
pable of influencing the shape of the nebula (the close and inter- 
mediate period binaries), a significantly smaller fraction than the 
observed fraction of non-spherical PNe which is ~80%. This indi- 
cates that mechanisms other than binarity are responsible for shap- 
ing PNe, such a s stellar rotation, magnetic fields or planets (see 
|Pe Marcoll2009l. for a review). 

Using our model and the observed number of red giant stars 
on the top one magnitude of the RGB, we predict that the num- 
ber of PNe in the central 25 square degrees of the LMC should 
be 548. This is in good a greement with the 541 ±89 PNe observed 
bv lReid & Parked ( hoOol) . This result suggests that nearly all low 
mass stars produce a PN in contrary to t he "Binary Hypothesis" 
dMoe & De Marcoll2006l ; lDe Marcoll2009h which suggests that bi- 
nary interaction is required to produce a PN. 

We have also predicted the orbital element distributions of se- 
quence E stars. The predicted light amplitude, orbital period and 
full velocity amplitude distributions agree well with the observa- 
tions, thus providing support for our modelling procedure. 
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5 CONCLUSIONS 

We have used the population of red giants that are ellipsoidal vari- 
ables in the LMC to predict the fates of binary red giants. In our 
standard model, we found 7-9% of PNe contain close binaries, 23- 
27% contain intermediate period binaries (P < 500 years) capable 
of influencing the shape of a PN, 46-55% contain wide binaries, 
5-6% are from merged stars and 3-19% are from single stars. The 
production rate of post-RGB stars, consisting of a He white dwarf 
and a low-mass secondary, is ~50% of the production rate of close 
binary PNNe. 

Our predicted fraction of close binary PNe (7-9%) is some- 
what lower than current observational estimates of 10-20%. Simi- 
larly, the observed orbital period distribution for close binary cen- 
tral stars of PNe has few stars with periods more than 1 day whereas 
the models predict that there should be a significant population of 
binaries with periods out to 100 days. The periods predicted by the 
models can be brought into better agreement with the currently ob- 
served periods by reducing a C E- The current observational samples 
are small and detection techniques work best for short period sys- 
tems. Larger samples and different detection techniques are needed 
to improve the significance of the current observational estimates. 

The mass distribution of all white dwarfs produced shows two 
peaks around 0.4 and 0.6 M corresponding to He and CO white 
dwarfs, respectively. The birth rate of He white dwarfs is predicted 
to be about 10% of the total white dwarf birth rate, in agreement 
with the observed birth rate ratio. 



REFERENCES 

Alcock C. et al., 1998, AJ, 115, 1921 

Bertelli G, Mateo M., Chiosi C.,Bressan A., 1992, ApJ, 388, 400 
Bertelli G, Girardi L., Marigo P., Nasi E„ 2008, A&A, 484, 815 
Blum R. D., et al., 2006, AJ, 132, 2034 
Bond H. E., 1994, ASPC, 56, 179 
Bond H. E., 2000, ASPC, 199, 1 15 
Bond H. E., Grauer A. D., 1987, fbs.xonf, 221 
Bond H. E., Livio M., 1990, ApJ, 355, 568 
Bond H. E., Ciardullo R., Meakes M. G, 1992, IAUS, 151, 517 
Bondi H., Hoyle F, 1944, MNRAS, 104, 273 
Bopp B. W, Stencel R. E., 1981, ApJ, 247, L131 
Bopp B. W., Rucinski S. M., 1981, IAUS, 93, 177 
Brown J. M., Kilic M., Brown W. R., Kenyon S. J., 2011, ApJ, 
730, 67 

Carlberg J. K., Majewski S. R., Patterson R. J., Bizyaev D., Smith 

V. V., Cunha K., 201 1, ApJ, 732, 39 
Chen X., Han Z., 2008, MNRAS, 387, 1416 
Cutri R. M., et al., 2003, yCat, 2246, 
Davis P. J., Kolb U., Knigge C. ., 201 1, ASPC, 447, 1 15 
de Kool M., 1990, ApJ, 358, 189 
de Kool M., 1992, A&A, 261, 188 
De Marco O., 2009, PASP, 121, 316 

De Marco O., Hillwig T. C, Smith A. J., 2008, AJ, 136, 323 

De Marco O., Passy J.-C, Moe M., Herwig F, Mac Low M.-M., 

PaxtonB., 2011, MNRAS, 411, 2277 
de Medeiros J. R„ Melo C. H. F, Mayor M., 1996, A&A, 309, 

465 



Predicting the fate of binary red giants 17 



de Ruyter S., van Winckel H., Dominik C, Waters L. B. F. M., 

Dejonghe H., 2005, A&A, 435, 161 
Dopita M. A., Meatheringham S. J., Webster B. L., Ford H. C, 

1988, ApJ, 327, 639 
Dorren J. D., Guinan E. F, 1984, NASCP, 2349, 391 
Driebe T., Schoenberner D., Bloecker T., Herwig F, 1998, A&A, 

339,123 

Duquennoy A., Mayor M., 1991, A&A, 248, 485 

Eggleton P. P., 1983, ApJ, 268, 368 

Frankowski A., Tylenda R., 2001, A&A, 367, 513 

Frogel J. A., Cohen J. G., Persson S. E., 1983, ApJ, 275, 773 

Frogel J. A., Mould J., Blanco V. M., 1990, ApJ, 352, 96 

Gielen C, et al., 2009, A&A, 508, 1391 

Girardi L., Bressan A., Bertelli G., Chiosi C, 2000, A&AS, 141, 
371 

Glebocki R., Stawikowski A., 2000, AcA, 50, 509 
Han Z., 1998, MNRAS, 296, 1019 

Han Z., Podsiadlowski P., Eggleton P. P., 1995a, MNRAS, 272, 
800 

Han Z., Eggleton P. P., Podsiadlowski P., Tout C. A., 1995b, MN- 
RAS,277,1443 

Houdashelt M. L., Bell R. A., Sweigart A. V., Wing R. F, 2000a, 
AJ, 119, 1424 

Houdashelt M. L., Bell R. A., Sweigart A. V., 2000b, AJ, 119, 
1448 

Hurley J. R., Tout C. A., Pols O. R., 2002, MNRAS, 329, 897 
Iben I., Jr., Renzini A., 1983, ARA&A, 21, 271 
Iben I., Jr., Livio M., 1993, PASP, 105, 1373 
Kamath D., Wood P. R., Soszyfiski I., Lebzelter T, 2010, MN- 
RAS, 408, 522 

Karakas A. I., Tout C. A., Lattanzio J. C, 2000, MNRAS, 316, 
689 

Keller S. C, Wood P. R., 2006, ApJ, 642, 834 
Kilic M., Stanek K. Z., Pinsonneault M. H., 2007, ApJ, 671, 761 
Kilic M., Brown W. R., McLeod B., 2010, ApJ, 708, 41 1 
Kiss L. L„ Bedding T. R., 2003, MNRAS, 343, L79 
Kouwenhoven M. B. N., Brown A. G. A., Portegies Zwart S. F., 

Kaper L., 2007, A&A, 474, 77 
Lebzelter T, Wood P. R., 2005, A&A, 441, 1 1 17 
Liebert J., Bergeron P., Holberg J. B., 2005, ApJS, 156, 47 
Lloyd Evans T, van Winckel H., 2007, IAUS, 240, 682 
Lutz J., Fraser O., McKeever J., Tugaga D., 2010, PASP, 122, 524 
Marsh T. R., Dhillon V. S., Duck S. R., 1995, MNRAS, 275, 828 
Massarotti A., Latham D. W., Stefanik R. P., Fogel J., 2008, AJ, 

135, 209 

Mastrodemos N., Morris M., 1999, ApJ, 523, 357 

Meixner M., et al., 2006, AJ, 132, 2268 

Mikolajewska J., 2007, BaltA, 16, 1 

Mikolajewska J„ 201 1, arXiv, |arXiv:1110.2361| 

Miller G. E., Scalo J. M., 1979, ApJS, 41, 513 

Miszalski B., Acker A., Moffat A. F. J., Parker Q. A., Udalski A., 

2009, A&A, 496,813 
Miszalski B., Corradi R. L. M., Jones D., Santander-Garcfa M., 

Rodriguez-Gil P., Rubio-Dfez M. M., 201 1, apn5.confE.328M 
Moe M., De Marco O., 2006, ApJ, 650, 916 
Morris S. L., 1985, ApJ, 295, 143 

Nicholls C. P., Wood P. R., Cioni M.-R. L., Soszyhski I., 2009, 

MNRAS, 399, 2063 
Nicholls C. P., Wood P. R., Cioni M.-R. L., 2010, MNRAS, 405, 

1770 

Nie J. D„ Zhang X. B., Jiang B. W„ 2010, AJ, 139, 1909 
Nikolaev S., Weinberg M. D., 2000, ApJ, 542, 804 



Paczyhski B., 1971, ARA&A, 9, 183 

Paczynski B., 1976, IAUS, 73, 75 

Politano M., Weiler K. P., 2007, ApJ, 665, 663 

Raghavan D., et al., 2010, ApJS, 190, 1 

Rebassa-Mansergas A., et al., 2008, MNRAS, 390, 1635 

Rebassa-Mansergas A., Nebot Gomez-Moran A., Schreiber 

M. R., Girven J., Giinsicke B. T, 2011, MNRAS, 413, 1121 
Refsdal S., Weigert A., 1970, A&A, 6, 426 
Reid I. N., Gizis J. E., Hawley S. L., 2002, AJ, 124, 2721 
Reid W. A., Parker Q. A., 2006, MNRAS, 373, 521 
Reimers D., 1975, MSRSL, 8, 369 
Rieke G. H., Lebofsky M. J., 1985, ApJ, 288, 618 
SalpeterE. E., 1955, ApJ, 121, 161 
SokerN., 1997, ApJS, 112, 487 
Soker N., Subag E., 2005, AJ, 130, 2717 
Soker N., Rappaport S., Haipaz A., 1998, ApJ, 496, 842 
Soszyfiski I., 2007, ApJ, 660, 1486 
Soszynski I., et al., 2004, AcA, 54, 347 
Stasinska G., Tylenda R., 1990, A&A, 240, 467 
Tout C. A., Eggleton P. P., 1988, MNRAS, 231, 823 
Tout C. A., Pols O. R., Eggleton P. P., Han Z., 1996, MNRAS, 

281, 257 

Udalski A., Szymahski M., Kubiak M., Pietrzynski G., Soszyhski 

I., Wozniak P., Zebruh K, 2000, AcA, 50, 307 
van Winckel H., 2007, BaltA, 16, 1 12 
van Winckel H., et al., 2009, A&A, 505, 1221 
Vassiliadis E., Wood P. R., 1994, ApJS, 92, 125 
Webbink R. F, 1976, ApJ, 209, 829 
Webbink R. F, 1984, ApJ, 277, 355 
Wood P. R., et al., 1999, IAUS, 191, 151 
Wood P. R., Faulkner D. J., 1986, ApJ, 307, 659 
Wood P. R., Olivier E. A., Kawaler S. D., 2004, ApJ, 604, 800 
Yungelson L. R., Tutukov A. V., Livio M., 1993, ApJ, 418, 794 
Zhang C. Y., Kwok S., 1993, ApJS, 88, 137 
Zijlstra A. A., 2007, BaltA, 16, 79 

Zorotovic M., Schreiber M. R., Gansicke B. T, Nebot Gomez- 
Moran A., 2010, A&A, 520, A86 



